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ABSTRACT 

■ The existence of flux-ratio anomalies between fold and cusp images in galaxy-scale strong- 

J> ' lens systems has led to an interpretation based on the presence of a high mass-fraction of 

•/^ , cold-dark-matter (CDM) substructures around galaxies, as predicted by numerical N-body 

simulations. These substructures can cause large perturbations of the image magnifications, 
leading to changes in the image flux ratios. The flux-ratio anomaly is particularly evident in 
the radio-loud quadruple gravitational lens system CLASS B2045-I-265. In this paper, new 
high-resolution radio, optical, and infrared imaging of B2045-I-265 is presented which sheds 
I more light on this anomaly and its possible causes. First, deep Very Long Baseline Array 

■~-«^ . observations show very compact images, possibly with a hint of a jet, but with no evidence 

' for differential scattering or scatter broadening. Hence, the flux-ratio anomaly is unlikely to be 

caused by refractive scattering in either the Milky Way or the lens galaxy. Second, optical and 
Q ■ infrared observations with the Hubble Space Telescope and through Adaptive-Optics imaging 

^ ' with the W. M. Keck Telescope, show a previously undiscovered object - interpreted as a 

, (tidally disrupted) dwarf satellite based on its colours and slight extension - between the 

main lens galaxy and the three anomalous flux-ratio images. Third, colour variations in the 
early-type lens galaxy indicate recent star-formation, possibly the result of secondary infall 
of gas-rich satellites. A population of such galaxies around the lens system could explain the 
p\ ' previously discovered strong [O ll] emission. However, spiral structure and/or normal star 

H I formation in the lens galaxy cannot be excluded. In light of these new data, we propose a lens 

model for the system, including the observed dwarf satellite, which reproduces all positional 
and flux-ratio constraints, without the need for additional CDM substructure. Although the 
model is peculiar in that the dwarf galaxy must be highly flattened, the model is very similar 
to recently proposed mass models based on high-order multipole expansions. 

Key words: gravitational lensing - quasars: individual: CLASS B2045-I-265 - cosmology: 
observations 



1 INTRODUCTION 

The standard model of cosmological structure formation predicts 
that galaxies form hierarchically from the mergers of smaller-mass 
progenitors, leading to increasingly more massive galaxies. Ac- 
cording to this model - as seen in numerical simulations - the 
baryonic stellar and gas components of galaxies should be sur- 
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rounded by an extended "smooth" dark-matter component, contain- 
ing many dwarf satellites. In cold-dark-matter (CDM) only simu- 
lations, this ensemble of compact satellites is gener i cally termed 
"CDM substructure" jNavarro. Frenk & Whit 3 1 19961: 1 Moore et al.l 
1998). However, if one compares the density of these satellites 
in numerical simulations with those found around, e.g., the Milky 
Way or M31, a severe discrepancy is apparent in the number of ob- 
served dwarf satellites. There ar e 1-2 orders of magnitude fewer 
dwarfs observed than predicted jKlvpin et alJll999riMoore et al.l 
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Il999 h. If this discrepancy is genuine, then either these satellites 
do not exist and the ACDM paradigm has a major problem, or they 
consist purely of dark matter having had their baryons stripped dur- 
ing the early phases of the galaxy-formation process. 

Strong gravitational lensing offers a powerful method to detect 
and quantify the presence of mass substructure, be it luminous or 
dark. This detection could be achieve d through the deta iled mod- 
elling of Einstein ring emission (e.g. lKooDmar^l2005h . Alterna- 
tively, substructure may also be detected through its effe ct on the 
flux-ratios between fold or cusp lensed imag es (Ma o & SchneideJ 
1 19981 : iKeeton, Gaudi & PettersI [iool I2OO5I ; iBradac et al.l l2004h . 
When the source distance to the fold or cusp is much smaller 
than the Einstein radius, the sum of the image brightnesses 
(assuming they are point images and not affected by variabil- 
ity) add to zero when the sign of the image parities are ac- 
counted for (i.e. negative-parity images are taken to have nega- 
tive brightness and magnification). Most lens systems with fold 
and cusp images appear not to obey these relationships be- 
tween image brightnesses, which in principle are independent of 
the global mass model. These flux-ratio anomalies have subse- 
quently been interpreted as the first direct evidence for the exis- 
tence of C DM substructure around lensing galaxies ('Bradac et al.l 
2002l:lchi ba 2002; Dalai & Kochanek 2002; Metcalf & Zhao 2002 



Keeton et al. 2003; iKochanek & Dalall 120041 ; iBradac et alj 120041 : 



Keeton et al.»2005i) . However, the mass fraction of this s ubstructure 



appea rs to be on the high side when compared to models jMao et aEI 
l2004h and the precise interpretation of these intriguing results re- 
mains open for debate. 

In this paper, we present new optical, near-infrared and radio 
data on the gravitational lens system CLASS B2045+265, which 
has one of the largest flux-ratio anomalies of all known lens sys- 
tems. These data, and in particular the discovery of a previously 
unknown optically-bright dwarf-galaxy between the main lensing 
galaxy and the three cusp images, sheds new light on whether the 
anomaly is caused by differential scattering, genuine dark-matter- 
only satellites or possibly by luminous dwarf satellites. In Section 
[2I we give a brief description of the previous observations and mod- 
els for B2045+265, and discuss the possible reasons for the flux 
ratio anomaly. Our new high-resolution radio data are presented 
in Section [3] The optical and infrared data sets, which include 
new adaptive optics imaging with the W. M. Keck Telescope and 
archival Hubble Space Telescope data, are presented in Section |4] 
The radio, optical, and infrared data are combined and used to con- 
strain a new mass model for B2045+265, which includes the pres- 
ence of the newly-discovered dwarf galaxy, in Section [5] Finally, 
in Section [6] we discuss our results and present our conclusions. 
Throughout this paper, we adopt an Q.m = 0.3, f^A = 0.7 flat Uni- 
verse, with a Hubble parameter of Hq — 100 h km s~^ Mpc^^. 



2 CLASS B2045+265 



lFassnachtet"ZI ( ll999h report the discovery of B2045+265, which 
is one of the most intriguing gravitational lens systems to be found 
during the course of the Cosmic Le ns All-Sky Survey (CLASS; 
iMvers et al.l2003l ; lBrowne et al.l2003.) . This system comprises four 
radio-loud lensed images in a long-axis cusp configuration; the 
three brightest lensed images (A, B and C) being almost co-linear 
in a north-south direction, and very close together on the sky 
(0.84 arcsec separation). The fourth and faintest lensed image, D, 
is 1.91 arcsec from the brightest lensed image. A, and is ~9 times 
fainter at 5 GHz. The system also has emission from a fifth radio 



component, E, which lies between D and the other lensed images. 
Multi-frequency radio mapping from 1.4 to 15 GHz with the Very 
Large Array (VLA) and Multi-Element Radio Linked Interferomet- 
ric Network (MERLIN) found the radio spectra of the four lensed 
images to be nearly identical (ail4 ~ —0.42 where Su oc v"). 
However, the radio spectrum and flux-density of component E was 
not consistent with a fifth lensed image. Therefore, E is believed 
to be emission from an AGN (Active Galactic Nucleus) within 
the lensing galaxy. Follow-up VLBA (Very Long Baseline Array) 
imaging at 5 GHz detected unresolved emission from lensed im- 
ages A, B and C. However, the short 35-min observation was not 
sensitive enough to detect image D or radio component E. 

High resolution infrared imaging of B2045+265 with the Near 
Infrared Camera/Multi-Object Spectrograph (NICMOS) on board 
the Hubble Space Telescope found a lensing galaxy coincident 
with the emission from radio component E and strong point source 
emission from the lensed images A, B and C with an underly- 
ing gravitational arc. The weakest lensed image was not detected 
in the infrared. Follow-up spectroscopic observations of the lens- 
ing galaxy showed strong [O II] emission and Balmer absorption 
lines which together are consistent with an Sa late-type galaxy at 
2 = 0.867 ( lFassnachtet^ll99 9!). A single, broad and isolated 
emission line was found in the spectrum of the background quasar, 
which has been identified as Mg II at a redshift of z = 1.28. How- 
ever, the uncomfortably low source redshift with respect to the lens 
requires a massive deflector, with a large dark matter component, 
to cause the 1.91 arcsec image splitting observed in B2045+265 
(see Fassnacht et al. 1999). Recently, Hamanaetal. (2005) mea- 
sured the stellar velocity dispersion of the lensing galaxy to be 
(j„ = 213±23 km s'^ This is significantly higher than the char- 
acteristic velocity dispersion for the late-type lensing galaxy pop- 
ulation obtained from the CLASS gravitational lensing statistics 
(cr„ = 117+ti; IChaeetalJl2002l) , leading to some doubt over 
the Sa-type classification of the lensing galaxy. Fu rthermore, the 
joint lensing and dynamics analysis carried out by iHamana et al.l 
concluded that the overall density profile, p(r), of the lens must 
be shallower than isothermal (/3 = 1.58 where p(r) oc r~^), 
which may be due in part to a possi ble nearby group of g alaxies 
not accounted for in the modelling ( IFassnacht et al.l [19991) . Con- 
versely, if the s ource redshift is higher than the 1.28 reported by 
IFassnacht et al.L then the lens could have a lower mass and hence a 
steeper density profile, consistent with the measured stellar velocity 
dispersion. 

Attempts to model the system with a smooth mass distribu- 
tion have successfully reproduced the positions of the lensed im- 
ages to within the astr ometric uncertainties , but failed to account 
for the flux ratios (e.g. IFassnacht et ai]|l999 ). For a smooth global 
mass model, B is expected to be the brightest of the three merg- 
ing lensed images; however, image A is significantly brighter than 
image B. The flux ratio anomaly of B2045-I-265 can be best quanti- 
fied using the dime nsionless Rcusp parameter tBradac et al.ii20o3 ; 
lKeetonetalj2003l) , 



Rc 



PA+ ^J.B + IJ.C 

+ |a*s| + Imc 



(1) 



where fiA, fJ,B and fic are the magnifications of merging images 
A, B and C, respectively. For a true cusp lens system with a glob- 
ally smooth mass distribution, Rcusp should tend to as the source 
position approaches the cusp, i.e., the flux density of B should 
approach the sum of the flux densities of A and C. In the cas e 
of B2045+265, Rcusp = 0.506±0.013 dKoopmans et"aDl2003ah . 
which is one of the most extreme examples of a cusp flux ratio 
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Table 1. The flux-densities, positions and angular sizes of the B2045+265 radio components, as measured with the VLBA at 5 GHz. The errors in the flux 
densities include a 5 per cent uncertainty in the absolute flux-density calibration of the VLBA. The flux density ratios are measured relative to image A 
(Al-fA2). The position of component Al is 20''47"'20.28593'' +26°44' 02.691 1". 
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-790.4±0.1 
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0.9 


56 
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Figure 1. The naturally weighted 5-GHz maps of the B2045+265 lensed images (A, B, C and D) and lensing galaxy (E) taken with the VLBA on 1997 October 
17. Lensed image A is resolved into a compact core (Al) with a faint extended jet to the south-west (A2). The remaining lensed images are compact, although 
there is evidence for a slight extension in image C. The weak and compact radio emission from component E is consistent with an AGN embedded within 
the lensing galaxy, Gl. The contours in each map are (—3, 3, 6, 12, 24, 48, . . .) X rTmap, the rms map noise given in Table[T] The restoring beam size is 
3.4x 1.9 mas^ at a PA of —5.5°. The grey-scales are in units of mJy beam"'^. North is up and east is left. 



anomaly known jKeeton et al.ll2003l : lMetca3l2005al) . Furthermore, 
the large value for Rcusp is mostly due to images B and C both 
being much fainter (by 25-50 per cent) than wh at is predicted from 
a simple mass model jPobler & Keetonll2006l) . This could be ac- 
counted for if there were several CDM substucture clumps within 
the lensing galaxy or along the line-of-sight to images B and C. 
However, this interpretation may be inconsistent with the level of 
substructure predicted from ACDM simulations. 

Alternatively, the flux ratio anomaly may be due to propaga- 
tion effects, such as differential scattering or free-free absorption. 
Both of these effects are believed to be responsible for changing the 
properties of radio-loud lensed images in othe r gravitational lens 
systems (e.g. [W inn, Rusin & Kochanekl f2004 isiggs et alj |2004 
iMittal. Porcas & Wucknitz ,20(371) . However, free-free absorption 
and refractive scattering are frequency dependent and the radio 



spectra of the four le nsed images of B2045 -I-265 are nearly identical 
from 1.4 to 15 GHz (Fassn acht et al .119991) . Therefore, the cause of 
the flux ratio anomaly in B2045-I-265 is almost certainly gravita- 
tional in origin. 



3 RADIO OBSERVATIONS 

In this section, we present new high resolution radio observations 
of B2045-I-265, obtained with the VLBA. 



3.1 VLBA 5-GHz observation 

The snapshot VLBA imaging presented bv lFassnacht et 

was only sensitive enough to detect the three brightest lensed im- 
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ages (A, B and C) of B2045+265 at 5 GHz. Therefore, with the 
aims of detecting the weakest lensed image (D) at mas-scales, in- 
vestigating the nature of the fifth radio component (E), and search- 
ing for extended structure in all of the lensed images, we undertook 
a deep observation of B2045-I-265 with the VLBA at 5 GHz. The 
observation was phase referenced and carried out over 10.7 h on 
1997 October 17, with ~6.5 h spent on B2045-I-265. The data were 
recorded in the left circular polarisation through 4 IPs, each with a 
bandwidth of 8-MHz. A bit rate of 128 Mb s"^ with 2-bit sampling 
was employed. Phase refere ncing and fringe findi ng were carried 
out with JVAS B2048-H312 dWilkinson et"aLlll998h and 3C 454.3, 
respectively. The data were correlated at Socorro, where the indi- 
vidual IPs were divided into 16x0.5 MHz channels and the visibil- 
ities were averaged over 1-s time intervals. 

The data editing and calibration was performed within AIPS 
(Astronomical Image Processing Software) using the VLBAPIPE 
pipeline script. No further frequency or time averaging was under- 
taken during this process. This produced an effective field of view 
for the resulting maps of ~5 arcsec from the phase tracking centre, 
easily large enough to contain the full extent of the lens system, 
which has a maximum lensed image separation of 1.91 arcsec. The 
image deconvolution was carried out using imagr. We first created 
a shallow wide-field map and identified the locations of the five 
radio components of B2045-I-265. A small 80x80 mas^ field was 
then placed around each component. A process of cleaning with 
two cycles of phase-only self calibration, using a solution interval 
of 10 min, was carried out. Natural weighting was used through- 
out. The resulting restoring beam size was 3.4x 1.9 mas'^ at a PA 
of —5.5°. The image rms in each field was ~55 ^Jy beam"^. Pi- 
nally, elliptical Gaussian model components were fitted to the radio 
components using jmfit. 
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Figure 2. Adaptive optics imaging at 2.2 fim with NIRC2 on the W. M. 
Keck-II Tele.scope. The.se new high resolution data have detected for the 
first time a faint galaxy (G2) which may be responsible for the flux ratio 
anomaly observed in B2045+265. North is up and east is left. 



flux density is too high to be a core lensed imageQ We again find 
this to be the case, with component E having a higher flux density 
than image D. Assuming that component E is associated with the 
lensing galaxy, we calculate the rest-frame 1.4 GHz luminosity of 
the component E radio emission to be 2.9 x 10^^ W Hz~^, 
which is consistent with a weak AGN. The flat radio spectrum 
{oi]^ ~ —0.4 between 5 and 15 GHz, where Sv oc v°') and com- 
pact nature of component E at mas-scales are both in agreement 
with the classification as a radio-loud AGN. The properties of the 
elliptical Gaussian fitted to component E are also presented in Ta- 
ble^ 



3.2 Results 

The deconvolved naturally-weighted maps of the four lensed im- 
ages are presented in Fig.[T]and show a mixture of compact and ex- 
tended emission consistent with gravitational lensing. Image A has 
been resolved into two subcomponents, forming a core-jet struc- 
ture. Subcomponent Al, the core, is compact and is the dominant 
source of the 5 GHz emission from image A. The faint extended 
nature of subcomponent A2 is consistent with a radio jet. Further- 
more, A2 is extended toward the south west, which is perpendicular 
to the axis connecting image A to the lensing galaxy. In contrast, 
image B is unresolved in our deep VLBA map. This is not surpris- 
ing since the surface brightness of the lensed images is conserved 
by gravitational lensing and image B is fainter than image A (the 
Sb/Sa flux density ratio is 0.61). The emission from image C is 
well represented by a single compact elliptical Gaussian. However, 
it is evident from Fig.[T]that there is a slight extension towards the 
south in image C, which is presumably due to faint extended emis- 
sion from the lensed jet. We have successfully detected the fourth 
lensed image of B2045-I-265 for the first time at mas-scales. Image 
D is the weakest lensed image and has a compact structure. The 
results of fitting Gaussian model components to the B2045-I-265 
lensed images are presented in Table [T] 

We have also detected the fifth radio component of 
B2045-F265 with the VLBA. The naturally weighted map is also 
shown in Fig.[T] We fi nd the component E r adio emission to be both 
weak and unresolved. iFassnacht et al.l ( ll999i) argue that component 
E is most probably emission from the lensing galaxy because the 



4 OPTICAL AND INFRARED OBSERVATIONS 

In this section, we present new adaptive optics imaging acquired 
with the W. M. Keck Telescope, and archival data taken with the 
Hubble Space Telescope. 

4.1 W. M. Keck Telescope adaptive-optics imaging 

High resolution ground based imaging of B2045-I-265 at 2.2 /im {K- 
band) was obtained on 2005 July 3 1 with the Near InfraRed Cam- 
era 2 (NIRC2; Matthews et al. in preparation) behind the adaptive- 
optics bench on the W. M. Keck-II Telescope. The NIRC2 narrow 
camera was used throughout, which provides a field of view of 
10 x 10 arcsec^ and a pixel-scale of 9.94 mas pixel" ^. The data 
were taken in twenty three 180-s exposures with a small dither be- 
tween each to facilitate good sky background subtraction during 
the reduction stage. A sodium laser guide star was used to correct 
for the atmospheric turbulence, and a V = 14. 1 -mag natural guide 
star, 33 arcsec from B2045-I-265, was used for fast guiding. 

The data were reduced within IRAl(3 using a double-pass re- 
duction algorithm. The first pass is used to map the positions of 

^ Recall that component E also has a flatter radio spectrum than the lensed 
images. 

^ IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by AURA, 
Inc., under cooperative agreement with the National Science Foundation. 
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Table 2. The relative positions and (Vega) magnitudes of the optical and infrared components of B2045+265. Note that radio component E is coincident 
with Gl. The relative positions have been taken from the 2.2 /im adaptive optics imaging. The magnitudes have not been corrected for galactic extinction. 
Non-dectections are quoted at the 3 a level. 
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the objects on the sky, which are then masked out during the sec- 
ond pass where the usual sky subtraction with temporally adjacent 
frames is carried out. The bad pixel and cosmic ray masks gener- 
ated during this process were then used to construct weight maps. 
The individual exposure s were combined usin g the weight maps 
and the drizzle package ( Fruchter & Hook 2002), which also cor- 
rected the data for the geometric distortions across the NIRC2 cam- 
era. We found that the quality of the corrected images varied greatly 
over the dataset we obtained on B2045+265, with the point-source 
full width at half maximum (FWHM) ranging from 67 to 93 mas. 
Therefore, we also drizzled subsets of the data together, which in- 
cluded the 6 and 12 best individual exposures i.e. those with the 
sharpest point spread function (psf). 

In Fig. [2] we show our adaptive optics assisted A'-band image 
constructed using all twenty- three exposures of B 2045-1-265 (psf 
FWHIVI of ^--^93 mas). First, it is clear that all of the radio com- 
ponents found in the VLBA maps have infrared counterparts (note 
that radio component E has been labelled Gl in Fig. |2l. Strong 
point-source emission from the four lensed images has been de- 
tected. Also, the infrared gravitational arc associated with A, B and 
C has once again been observed. The lensing galaxy Gl, which is 
coincident with radio component E to within the positional uncer- 
tainties, has a circular morphology (axial ratio of b/a = 0.94±0.01 
at a position angle of —60.4° ±4.8° east of north) with no evi- 
dence of spiral arms. Further more, the surface b rightness profile of 
Gl, measured using GALFIT jPeng et al]|2002h . is consistent with 
a de Vaucouleurs model with an effective radius of 0.465 arcsec 
(= 1.65 h^^ kpc). This result is not too surprising since these K- 
band observations correspond to 1.2 /im in the rest frame of Gl, and 
are therefore sensitive to the older (red) stellar population contained 
within the bulge. Finally, and most importantly, we have also de- 
tected infrared emission from a previously unknown object (which 
we label G2), positioned between the main lensing galaxy Gl and 
the lensed images A, B and C. G2 is significantly fainter than Gl 
(the Gl /G2 flux ratio is ~ 64). In Section[5] we will test whether 
G2 is responsible for the large flux ratio anomaly between the three 
merging images of B 2045 -1-265. 

In Table|2]we report the relative positions and apparent magni- 
tudes for the infrared components of B2045-I-265. No photometric 
standard stars were observed during our adaptive optics observa- 
tions. Therefore, we compared the brightnesses of stars in our adap- 
tive optic s imaging and the calib rated NIRC AT-band imaging pre- 
sented bv lFassnacht et al.l ( 1 1999 1) to calculate the NIRC2 zeropoint 
magnitude - note that there were only two isolated stars common to 
both images. The absolute photometric accuracy resulting from this 
process is estimated to be no better than ~ 0. 15 mags. Also, the rel- 
ative photometric accuracy from our adaptive optics imaging will 



be limited by changes in the psf, and therefore the Strehl ratio, over 
time and position on the detector (e.g. Esslinger & Edmunds 199^. 
For Gl we state the total magnitude from the fitted de Vaucouleurs 
profile. A circular aperturqj was used to measure the magnitudes 
of the lensed images and G2, from a Gl subtracted image. 

4.2 Hubble Space Telescope imaging 

To investigate the nature of G2 further, we have retrieved high reso- 
lution optical and infrared imaging of B2045-I-265 from the Hubble 
Space Telescope archive. The optical data were taken with the Wide 
Field Planetary Camera (WFPC2) on 1998 December 04 and con- 
sist of two 500-s exposures through both the F555W and F814W 
filters (GO-6629; PI: Jackson). The lens system was positioned in 
the centre of the Planetary Camera (PCI), which has a field of view 
of 37 X 37 arcsec and a pixel scale of 46 mas pixel" \B2045-l-265 
was re-observed with NICMOS on 2003 September 6 (GO-9744; 
PI: Kochanek). Th ese new observations w ere much deeper than 
those presented by iFassnacht et al] jl999l) . and were designed to 
detect the fourth lensed image and provide high signal to noise ra- 
tio imaging of the lensing galaxy. The data were acquired with the 
FI60W filter and the NIC2 camera, which has a field of view of 
19 X 19 arcsec'^ and a pixel scale of 75 mas pixeP^. The imag- 
ing was taken in five 640-s and three 704-s exposures (total time 

5.3 ks). 

The WFPC2 and NICMOS datasets were downloaded from 
the Multimission Archive at Space Telescope and processed within 
IRAF using standard STSDAS pr ocedure s. In particular, the mul- 
tidrizzle package dKoekemoer et al.l2002h was used to combine the 
individual exposures in each filter, remove cosmic ray events and 
correct for the WFPC2 and NICMOS geometric distortions. Our 
final optical and infrared images of B2045-I-265 are presented in 
Fig. [3] Note that the WFPC2 imaging was not sensitive enough 
to detect all of the components of B2045-I-265. Image D was not 
detected at either F555W or F8I4W. Furthermore, G2 was not 
found in the F555W image and the main lensing galaxy Gl was 
only marginally detected. These non-detections are due in part to 
the short exposure times through both the F555W and F8I4W fil- 
ters with WFPC2, but are also due to the high extinction along 
the line-of-sight to B2045-I-26 5 owing to its low g alactic latitude 
{b = - \QA°;Av = 0.8 mags: ISchlegel et al.ll998b . The positions 
of the missing components in the F555W and F814W images have 
been marked by the dashed circles in Fig|3] 

^ The radii of all circular apertures used throughout this paper are set such 
that, for a point source, the measured magnitudes from aperture photometry 
and psf fitting would be identical. 
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Figure 3. Optical and infrared imaging of B2045+265 with WFPC2 and NICMOS on tlie Hubble Space Telescope. Tlie WFPC2 F555W and F814W data were 
not sensitive enough to detect all the components of B2045+265. A dashed circled has been placed over the positions of the missing components. North is up 
and east is left in each image. 



The structure of the lensing galaxy is significantly different in 
the F814W and F160W images. At F814W, there is clear evidence 
of an irregular morphology which would be consistent with emis- 
sion from either spiral arms or additional infalling dwarf galaxies. 
Alternatively, the observed structure could be artifacts of the data 
reduction process (we cannot rule out the possibility that these fea- 
tures are residual cosmic rays or bad pixels as only two undithered 
observations of B 2045 -1-265 were taken with WFPC2). Conversely, 
the lensing galaxy has a smooth surface brightness distribution at 
F160W, which is well fitted with a de Vaucouleurs profile. The fit- 
ted effective radius (0.411 arcsec) and axial ratio (0.92±0.01 at 
a position angle of —66.0° ±4.3° east of north) of Gl at F160W 
are both consistent with the results found using the A'-band adap- 
tive optics imaging presented in Section 143] The lensed images 
of the background quasar are unresolved point sources at F555W, 
F814W and F160W (note that the gravitational arc associated with 
A, B and C was found at F160W after subtracting the lensed im- 
ages - but is not shown here). The flux ratio anomaly between the 
three merging lensed images is also present in the optical and in- 
frared (e.g. at F160W the flux-ratios are Fb /Fa = 0.65±0.04 and 
Fc I Fa ~ 0.89±0.04), which further adds to the case for the amo- 
maly being gravitational in origin. In an attempt to characterise the 
surface brightness properties of G2, we have fitted the F160W data 
for G2 with a simple Gaussian profile. We find the G2 emission to 
be slightly extended and elongated toward the lensed images. The 
fitted Gaussian for G2 has an axial ratio of 0.56±0.18 at a posi- 
tion angle of 106.3°±21.1° east of north. The G1/G2 flux ratio at 
F160W is ~ 67. 



In Table[2]we give the optical and infrared magnitudes of GI, 
G2 and the lensed images. For GI, we report the total magnitude for 
the best fitting de Vaucouleurs profile (for the F555W and F814W 
datasets, we have held fixed the parameters of the de Vaucouleurs 
profile fitted to the F160W data, and allowed only the total magni- 
tude to vary). Point-sources were used to determine the magnitudes 
of the lensed images through each filter. For G2, we have measured 
the magnitude within a circular aperture from an image with the 
lensing galaxy and lensed images subtracted. 



4.3 The nature of G2 

We have labelled the newly detected object G2 which implies that it 
is a second galaxy along the line-of-sight to B2045-I-265. However, 
given the low galactic latitude of B2045-I-265 there is the possi- 
bility that G2 may be a star. The strongest evidence for G2 being 
a galaxy comes from the extension and elongation detected in the 
FI60W imaging. We have compared the size of G2 with the size of 
10 stellar objects in the FI60W field (the F160W imaging has the 
highest signal-to-noise ratio for G2 of all the available datasets); 
the average deconvolved FWHM of the stars is 0.73 pixels (rms of 
0.22), whereas for G2 the deconvolved FWHM is 1. 71 ±0.30 pix- 
els. 

Measuring the spectrum and obtaining a redshift for G2 would 
provide conclusive evidence that it is inde ed a galaxy. Ex t ensive 
spectroscopic observations of the system bv lFassnacht et al.l j 1 9991) 
were taken at a position angle (112°) which also placed the slit 
over G2. The resulting spectrum showed strong [O II] emission at 
the same redshift as lensing galaxy GI, which could be due to star- 
formation associated with either GI or G2. However, these ground 
based spectroscopic observations did not spatially resolve G2 from 
GI, and therefore the exact source of the [O II] emission is un- 
known. 

In an attempt to constrain the redshift of G2, we have cal- 
culated a photometric re dshift using the spectral energy di stribu- 
tion fitting code HYPERZ ( Bolzonella. Miralles & PeIIol200oh . This 
code compares the observed spectral energy distribution with an 
ensemble of theoretical/observed template spectra to establish the 
most probable redshift and spectral type. The optical and infrared 
data (C) presented in Table (2] after corr ecting for galactic ex tinc- 
tion using the Eb~v reddening maps of ISchlegel et alj (Il998h and 
the J? V -dependent extinction law of ICardelli, Clayton & Mathij 
( Il989l) . provided the observed photometric properties of G2. The 
template spectra (T) for elliptical (instantaneous burst, E/SO), spi- 
ral (Sa-Sd) and irregular (Ir) galaxy ty pes with 51 different stel - 
lar population ages were taken from the lBruzual & CharlotI jl993h 
GISSEL98 spectral library. The resulting G2 probability-redshift 
distribution for each template, p(C|2,T), is given in Fig.|4] 

There are two clear peaks in the probability-redshift distribu- 
tion for each galaxy template; the first is at Zphot ~ 0.86 and the 
second, which is much broader in redshift space, is at Zphot ~ 4.50. 
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This double-peaked probability-redshift distribution is a typical re- 
sult from photometric redshift codes when there is insufficient cov- 
erage of the observed spectral energy distribution. This is because 
strong spectral features such as the 4000 A and Lyman breaks are 
used to determine the photometric redshift and it is not always clear 
which break has been detected from only a few photometric mea- 
surements. To overcome this, we have also applied a Bayesian prior 
using the apparent magnitude (mp) of G2. Here, we follow the pro- 
cedure described bv lBem'tezl (l200Q^ where redshift distributions per 
spectral type are calculated from Hubble Deep Field-North data 
and used to establish the relative redshift-probability for a given 
template and apparent magnitude. This process essentially places a 
luminosity limit per spectral type e.g. there are fewer L* ellipticals 
at redshift 5 than at redshift 1, thus downweighting any high red- 
shift peak. In Fig. [5] we show the probability-redshift distribution 
for G2 using the prior on the apparent magnitude, p(z;|C,mo). We 
see that for the elliptical and spiral spectral types the high redshift 
peak is completely removed. For the irregular galaxy template the 
high redshift peak is significantly lowered and is now at a similar 
level to the low redshift peak. This analysis demonstrates that G2 
is unlikely to be a field galaxy at z = 4.50 because the implied lu- 
minosity is too large (Mt — —21.95; 4 L*). Furthermore, if G2 
were at z = 4.50 then we would expect it to also be gravitationally 
lensed by Gl. However, there is no evidence of a counter lensed 
image in either the Keck adaptive optics or the Hubble Space Tele- 
scope imaging. 

The best fitting template for G2 is an elliptical galaxy (instan- 
taneous burst) at Zphot ~ 0.i6t.Q'2i (68 per cent confidence lim- 
its) with a young 0.5 Gyr stellar population and low dust content 
(Av = 0). However, some caution is expressed over the proper- 
ties of the best fitting template since there are known to be strong 
degeneracies between stellar age and dust extinction. It is intrigu- 
ing that the most probable photometric redshift corresponds almost 
exactly to the redshift of the lensing galaxy Gl (z = 0.867). How- 
ever, this may just be a fortunate coincidence. For example, the 
derived photometric redshift for Gl is 0.46lo'26 (68 per cent con- 
fidence limits), which is only just consistent with the spectroscopic 
redshift. Therefore, further photometric measurements of G2 will 
need to be carried out to establish the full observed spectral energy 
distribution, and produce a more robust estimate of the photometric 
redshift. 




Figure 4. The G2 photometiic redshift probability distribution for elliptical, 
spiral and irregular galaxy templates found using HYPERZ and the optical 
and infrared data presented here. Due to the limited number of photometric 
data points available, there is a two peaked probability-redshift distribution. 
The photometric redshift of the two peaks are Zphot ~ 0.86 and 4.50. 




Figure 5. The photometric redshift probability distribution for G2, with the 
addition of a Bayesian prior on the apparent magnitude. The best fitting 
template for G2 is an elliptical at Zphot = 0-86^g'2^ ■ 



5 GRAVITATIONAL LENS MODEL 

In this section, we test whether G2 could be responsible for the flux- 
ratio anomaly observed between the three merging cusp images by 
constructing a new mass model for B2045-I-265. We assume that 
G2 is a dwarf satellite, or its remnant, at the redshift of Gl, which 
is consistent with the results from the spectral energy distribution 
fitting above. 

As galaxy Gl is well-fitted by a de Vaucou leurs profile, we as- 
sume it to be an early-type elliptica l galaxy. In lTreu & Koopmanj 
( |2004|) and iKoopmans et alj jBo^ it was shown, based on com- 
bined gravitational lensing and stellar kinematic data, that an 
isothermal mass model (i.e. p{r) cx r~^) often provides an ex- 
cellent description of the smooth density profiles of massive early- 
type galaxies out to a redshift of one, within their inner 0.3- 
5 effective radii. Even th ough the stell a r velo city dispersion of 
213±23 km s~^ found bv lHamana et al.l ( |2005[) indicates a some- 
what more shallo w density profile than i sothermal, the source red- 
shift identified bv lFassnacht et al.l ( 1 19991) is still disputed. A higher 



source redshift would lead to a lower lensing galaxy mass and 
consequently a steeper density profile in order to match the ob- 
served stellar velocity dispersion. B20 45+265 also exh ibits a possi- 
ble overdensity of nearby galaxies (e.g. lFassnacht et aL 1999) about 
~10 arcsec west of the lens system. 

We model the mass distribution of the system with two singu- 
lar isothermal ellipsoids (SIE), one for Gl and one for G2, plus an 
external shear. The constraints on the mass model are provided by 
the positions of the lensed images and the radio-loud AGN within 
Gl (0.1 mas accuracy) obtained from the new VLB A observations 
(see Table [T), supplemented by the average flux-ratios obtained 
from a M ERLIN Key-Proeram m onitoring program (~2 per cent 
accuracy; ' Koopmans et alj|2003ah . In order not to underconstrain 
the model, we fix the mass centroid of Gl and G2 to the brightness 
centroid measured from the Keck adaptive optics imaging (see Ta- 
bleO. 

The SIE mass models of Gl and G2 are free to change their 
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Table 3. Lens model parameters for B2045+265. Given are the lens strength 
(6), position angle (0) and axial ratio (g), as well as the external shear 
strength and its position angle (7, 9)ext- The source position and normal- 
ized flux are given in the last three rows (the source flux 5s is normalized, 
assuming Fp^ = 1.). The second column indicates the most likely value, 
whereas the third column indicates the (marginalized) median value from 
the MCMC simulation and the 68 percent confidence levels. 
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0.66±0.01 
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0.215±0.002 




6cxt 
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-0.641 


-0.641±0.001 


arcsec 




-0.654 


-0.654±0.001 


arcsec 


Fs 


0.104 


0.104±0.003 





axial ratios, position angles and lens strengths. The shear strength 
and position angle are also free to change. The total number of 
free parameters is 11, whereas the number of constraints is 12. 
The number of degrees of freedom is therefore NDF = 1 and 
the models should result in ~ 1. The model parameters are 
quickly found by optimizing the pena lty function, using the 
code described in iKoopmans et all ( l2003bl) . resulting in a best-fit 
model with xVNDF = 1.90 The resulting model of the system is 
shown in Fig.|6]and the best fitting parameters are listed in Table[3] 
with non-linear 68 per cent confidence levels obtained through a 
Markov-Chain Monte-Carlo simulation of 10^ samples. Note that 
the confidence limits quoted in Table[3]are the formal uncertainties 
and the systematics associated with the choice of model are much 
larger. 

We note several points. First, the model requires a significant 
external shear of ^--^0.2. Although large, it is not unexpected given 
the overdensity of nearby galaxies in the western direction of the 
shear position angle. However, the position angle of Gl is offset 
by '-^90 degrees from that of the external shear, which could point 
to a compensating effect, i.e. the lens potential is more boxy than 
allowed by the elliptical mass model. Hence, part of the external 
shear might be the result of a non-elliptical halo around Gl. Given 
the very tight constraints on the image positions and flux-ratios, 
these results are very robust (under the assumption that the mass 
model is chosen correctly). Second, the mass axial ratio of G2 is 
very small (~0. 13). Such values can only be expected from either 
edge-on disk system s or possibly a tidally disrupted dwarf galaxy 
(see e.g. Li et al. 2006). Whether this is the correct interpretation of 
the model results remains to be tested with additional, deeper, ob- 
servations. We note that although the mass axial ratio is more elon- 
gated than the light axial ratio measured from the F160W imag- 
ing presented in Section |4~2l the position angle of the G2 elonga- 
tion is consistent with the position angle of the light distribution 
(observed qG2 = 0.56±0.18 and 602 = 106.3°±21.r). Third, 
the resulting critical curve (Fig. [6} is similar to that produced by 



* The density slope found bv iHamana et alj i2005h . f3 = 1.58, gives 
= 5.3. Leaving the density slope as a free parameter results in /3 = 1.8 
with = 0.0 (this case has NDF = 0). 
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Figure 6. The B2045+265 mass model which reproduces the relative posi- 
tions and flux density ratios for the lensed images (A, B. C and D) without 
the need for additional cold dark matter substructure. The mass model con- 
sists of two singular isothermal ellipsoids representing Gl and G2, with an 
external shear The background source position is marked by the star. The 
lens critical curve is shown by the thick (red) line. The thin (blue) line shows 
the source plane caustics. The model parameters are given in Table[3] 



ICongdon & KeetonI (2005; their Figure 4) near images A-C. Even 
though they use a completely different method to model the system, 
the protrusion of the critical curve near the cusp images is similar to 
what a dwarf disk galaxy, or a tidally disrupted dwarf galaxy (e.g. 
iLi et alj|2006h . would cause when situated near G2. 



6 DISCUSSION & CONCLUSIONS 

We have presented high-resolution radio, optical and infrared im- 
ages of the gravitational-lens system CLASS B2045-I-265, obtained 
with the VLBA, Hubble Space Telescope and Keck (adaptive- 
optics) Telescope, respectively. In particular, the data were an- 
alyze d to examine the nature of the strong flux-ratio anomaly 
(e.g. iMao & Schneidej 1 19981) between the three cusp images. 
This anomaly has been ascribed to the prese nce of CDM sub- 
struct ure around the main lensing galaxy (e.g. Dalai & Kochanekl 
2002h. as predicte d from numerical simulations jMoore et al 1 19981 : 
Klvpin et aljl999 ). As such, these anomalies could be an important 
probe of the CDM paradigm. 

We have shown the following. First, the lensed radio images, 
because of their compactness in the VLBA imag es, are not strongl y 
affected by differential scattering at 5 GHz (e.g. lBiggs et alj2004t) . 
even though t he images are sc intillating due to the Galactic ion- 
ized medium jKoopmans et al.i r2003a). In addition, the flux-ratio 
anomaly observed at radio wavelengths (around 5 GHz) persists 
over at least eight months ( Koopmans et alji2003a) . The flux-ratio 
anomaly is therefore most likely due to a genuine perturbation of 
the smooth mass model of the lensing galaxy. Second, similar flux- 
ratios between the radio, optical and infrared further support that a 
mass perturbation in the lensing galaxy is causing the anomalies in 
this system, each of which individually could be affected by differ- 
ent processes (e.g. dust and microlensing in the optical, scattering 
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in the radio, etc.) but unlikely in a way that is similar in magnitude 
over this wide range of wavelengths. 

This leaves open the question what the true nature of the mass 
perturbation is. Are they caused by purely dark CDM satellites (e.g. 
[Oalal & Kochanek 2002; Kochanek & Dalai 2004), or by luminous 
CDM substructure in the form of normal dwarf satellites as we see 
them around the Milky Way and other galaxies? 

In the present Universe, luminous dwarf satellites might not 
be su fficient to explain these anomalies (e.g. Dalai & Kochanek 



1 20021) . However, we expect their numbers to be more abundant 
at higher (lens) redshifts, simply because many of them will have 
merged with their hosts between the lens r edshifts and z ^ 0. For 
example, the anomaly in MG 2016+112 jKoopmans et al.l [20021) 
is most likely caused by a nearby dwarf satel lite detected through 
its [O I I] emission at the lens redshift (Koopmans &Treull2002l; 



see al so lKochanek. Schneider & Wambsganssi2004l and lChen et al.l 
l2007h . Another example is MG 0414+053 where the addition 
of a lu minous dwarf galaxy fou nd close to the main lensing 
galaxy jSchechter & Moorg Il993h results in a model which re- 
produce s the observed positions and flux- densities of the lensed 
images jTrotter et alj|2000l : IRos et alj[200ol) . However, the case of 
B2045+265 i s one where the flux-ratio anomaly is very strong (e.g. 
iKeeton et al. |[2003; Metcalf 2005a) and therefore of particular in- 
terest when investigating the nature of these anomalies. 

As in the cases of MG 2016+1 12 and MG 0414+053, we find 
evidence for the presence of at least one dwarf satellite G2 near 
the main lensing galaxy Gl (Figs. [2] and [3j in the Hubble Space 
Telescope images with WFPC2-F814W and NICMOS-F160W, and 
in the ^-band Keck adaptive optics images. G2 resides in between 
the lensing galaxy and the cusp images. In addition, the lensing 
galaxy shows significant structure in the F814W-F160W colour im- 
ages, suggesting the presence of starforming dwarf satellites. The 
[G II] emission detected from this system supports this interpreta- 
tion, but is not conclusive because it cannot be spatially matched (as 
in MG 2016+112). If the half-dozen distinct peaks in the F814W 
image around Gl are all dwarf satellites and not image artifacts, 
then this number is somewhat larger than the expected number of 
about three observable dwarf satellites, given that G2 is ~ 1.5 mag 
above the 3-a limit for detectio n and the local luminosity function 
of dwarf satellites is oc L~^^ jSmith. Martinez & Grahanj|2004l : 
we normalise the integral over the luminosity function above the 
luminosity of G2 to unity). Having built the case that G2 can be 
a luminous dwarf satellite of the lensing galaxy Gl, we have pre- 
sented a new lens model for B2045+265, taking this second bright- 
est group member into account (G2). The result indicates that the 
satellite mass distribution must be significantly flattened, which 
could either be due to an edge-on disk or possibly a tidally dis- 
rupted dwarf galaxy being accreted on to Gl (see e.g. lLi et al.l2006l 
for the effect of these streams on lens image flux-ratios). We realize 
that this model is speculative at present. However, the lens-model 
critical curves are remarkably simi lar to an independen t multi -pole 
expansion model of this system bv lCongdon & Keetonl ( |2005|) . 

It now seems likely that three out of the seven gravitational 
lens systems used to statistically detect dark matter substructure 
around lens galaxies jDalal & Kochane3l2002h show evidence for 
CDM substructure in the form of luminous dwarf satellites. Fur- 
thermore, this luminous substructure appears to be more mas- 
sive than expected when compared with the results from numer- 
ical simulations. For example, in the cases of B2045+265 and 
MG 2016+112 the luminous dwarf satellite constitutes about 1 per 
cent of the total halo mass, and for MG 0414+053 about 0.3 per 
cent. This mass-fraction, which is sensitive to the region contained 



within the lensed image positions, seems somewhat larger than 
to what is predicted for the inner few kpc of ga laxy-scale haloes 
from numerical simulations i.e. ^ 0.5 per cent jMao et aklliooi; 
iDiemandT Kuhlen&Madaull2007h . This discrepancy between the 
observed and simulated mass-fraction could be due to an over- 
merging problem at the cores of simulated haloes or an over- 
estimation of the effect tidal stripping by baryons has on dark mat- 
ter substructures in the numerical simulations. Alternatively, there 
could be additional unseen CDM substructure contributing to the 
flux-ratio anomaly which is not associated with the main lensing 
halo, but lies along the line of sig ht to the lensed images of the 
background quasar ( lMetcalj|2005br) . 

In summary, the early-type lens galaxy in the system 
B2045+265 shows evidence for additional mass structure in the 
form of luminous dwarf satellites. The system also shows one of 
the most anomalous flux-ratios between its three cusp images of 
any lens system. Although we have not definitely shown that the 
dwarf satellite (G2) is responsible for the anomaly, the presence 
of significant structure in the brightness of Gl and the presence 
of G2, makes it more likely than previously assumed. A simi- 
lar explanation was found for the anomalies i n MG 2016+112 
jKochanek et"ai]|2004l : 'Koopmans & Treu"2002'; 'Chen et alj |2"007h 
and MG 0414+053 feotter et al. 2000; Ros et al. 2000). Hence, 
deep high-resolution optical and infrared observations of each 
anomalous system is highly desirable to settle the question whether 
any or all anomalies are due to either luminous dwarfs or purely 
dark CDM substructures. 
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